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Abstract 

The problem of choosing an optimal strategy for moving in the rain has attracted 
considerable attention among physicists and other scientists. Taking a novel 
approach, this paper shows, by studying simple shaped bodies, that the answer 
depends on the shape and orientation of the moving body and on wind direction 
and intensity. For different body shapes, the best strategy may be different: in 
some cases, it is best to run as fast as possible, while in some others there is an 
optimal speed. 

(Some figures may appear in colour only in the online journal) 



1. Introduction 

Is it better to walk or run in the rain? Almost everyone has at least once faced this question; many 
discussions on this topic can be found on the web, with various conclusions. The problem has 
also been addressed by some physicists, mathematicians, engineers and meteorologists. Some 
studies [1] consider vertically falling rain. Others [2, 3] take into account the possibility of wind 
with the same direction of motion. Others still [4-8] also consider a cross (i.e. perpendicular 
to the path) component of the wind. Since the human body has a very complex shape, a 
much simpler form, usually a parallelepiped, has been considered, with the implicit or explicit 
assumption that shape is not a crucial factor, and that results obtained for a parallelepiped can 
be quite easily generalized to bodies of any shape. 

Qualitatively speaking, the results found up to now can be summarized as follows. 

• For rain falling vertically, the best strategy is to run as quickly as possible. The same is 
also true for motion into the wind. 

• For motion downwind, there may be an optimal speed, which equals the component along 
the direction of motion of the wind velocity. This happens only if the ratio between the 
cross-section of the body perpendicular to the motion and the horizontal one is large 
enough; otherwise, the best choice is again to run at the maximum speed one can reach. 

In essence, in the literature so far published on this subject there are two main conclusions: 
for the existence of an optimal speed the wind must come from behind and the optimal speed, 
when it exists, equals the component of the wind velocity along the path. In this paper, we 
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shall show that the answer to the question actually depends on many factors, mainly on the 
shape of the body and its orientation. Even though a tailwind is a favourable condition for the 
existence of an optimal speed, in some cases we can have one even with a headwind, and in 
general its value is not equal to that found in previous studies. 

The paper is addressed mainly to undergraduate students and their teachers, as well as to 
anyone intrigued by the problem. The subject requires good geometrical visualization, a good 
command of vector calculus and involves many concepts from various fields; as such, it could 
be useful in physics courses at university level. 

2. Definition of the problem 

The question stated at the very beginning of this paper is too vague, and should be better 
defined. 

Let us consider a fixed track that has to be covered in the rain. We shall assume the 
following. 

(1) The ground is horizontal. 

(2) The path is rectilinear. 

(3) The rain is uniform, at least along the path. 

(4) The rain is constant, at least for the time needed to cover the path. 

(5) The body motion is rigid and translatory. 

This last assumption is clearly a quite crude approximation of a human body, but it is 
necessary to make the problem manageable. The conditions (2), (3) and (4) are not as restrictive 
as they may look at first sight: if they are not satisfied on the whole path, we can always split 
it into shorter parts. 

3. Introduction of conceptual tools 

Let us suppose that we have to go straight from point A to point B. We choose a Cartesian 
reference system with the x axis oriented from A to B, and the z axis vertical. Let u be our 
velocity; in this frame, u = (u, 0, 0), with u > 0. The problem is to find the value of u, if any, 
that will minimize the mass of water impinging on our body along the path. 

Let v = (v x , v y , v z ) be the rain velocity. The vertical component, v z , depends on the drop 
size. We shall assume that drops fall at their terminal speed, so that their horizontal velocity 
equals that of the wind (we do not restrict ourselves to a wind along the x direction, but we 
assume a horizontal wind). The sign of v y simply tells us if we get wet on our left or right side. 
When v x is positive, we shall speak of a tailwind, otherwise of a headwind, irrespective of the 
presence of a cross component v y . 

The plane identified by u and v plays an important role in the problem: we shall call it ic ; 
an axis perpendicular to this plane will be denoted by f . 

Let us call p the ratio between the mass of water drops that are found, at some instant, 
within a given volume, and the volume itself (note that p is not the water density). We can 
now define a vector j that we shall call rain density: 

jo = pv. (1) 

It is easy to recognize the analogy between p and charge density, and between jo and the 
current density vector in electromagnetism. Our approach will freely exploit conceptual tools 
taken from this field. Assumptions (3) and (4) imply p and v — and consequently jo — to be 
uniform and constant. 
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Figure 1 . The relative orientation of v, u, j and the plane n, which contains all of them. 



When moving, we shall perceive an apparent rain density j, which will differ from jo: 

jO) = p(v - u) = p(v x - u, v y , v z ). (2) 

The j vector changes, in intensity and direction, with u. Its direction rotates, when u 
increases, in the plane n, drifting away from the positive x direction; we can easily see this 
effect on train or car windows. In astronomy, this effect is known as stellar aberration. Sailors 
know well that the direction from which the apparent wind comes approaches the bow as boat 
speed increases. We can state this concept more precisely by saying that the angle 6 between 
j and u increases with w, approaching 180°. It is easy to show that 



v* + v\ 

tan<9 = 1 . (3) 

v x -u 

We can note that jo is the value of j corresponding to u = 0; we shall call it the initial 
value of j. Similarly, the angle between jo and u will be denoted by 0q. 
Figure 1 shows the relative orientation of v, u, j, and the plane it. 

Using once more the analogy with electromagnetism, it is natural at this point to introduce 
the concept of rain flux. If we adopted without any change the flux definition for a vector field, 
we should define rain flux as the surface integral of the j field: 

d>( M ) = j£j-dA, (4) 

where the integration is over the whole surface of the body. But we know that, in the absence 
of field sources or wells, the integral over a closed surface vanishes: thus the integral (4) must 
be modified for our purposes. 

A first change is to consider the modulus of the dot product j • dA. But this is not 
satisfactory, as every field line contributes twice to the flux: when it enters the body and when 
it leaves. So we shall restrict the integration to what we shall call the 'wet surface' of our body, 
S w , that is the surface where field lines enter the body. If the body has a complex shape, so that 
there are some field lines that enter, go out and enter again, the wet surface obviously includes 
only those parts where field lines enter the body for the first time. We then define rain flux as: 

*(*0 = jf IJ' dA l- (5) 

Note that, while S is a closed surface, S w is always an open one, as the different notation 
of the integral remind us. We shall assume that it maintains the same orientation as S (which, 
by convention, is oriented from inside to outside). 
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Rain flux is the ratio between the mass of water impinging on the body during a given 
time interval and the time interval; its SI measure unit will be kg s _1 . It depends on the body 
speed for two reasons, because j changes with u and because, due to the rotation of j, the 
surface of integration may change. 



4. Formalisation and general discussion of the problem 

We can now, at least in theory, calculate the water mass, m, hitting the body during motion in 
the rain. The time interval, At, needed to travel a distance L at a constant speed u is At = L/u, 
so that m is given by 

m(u) = O At = - { |j -dA|. (6) 

u Js w 

While this equation is adequate to study the behaviour of solid bodies enclosed in plane 
surfaces (e.g. a parallelepiped), in the presence of curved surfaces it is in general more 
convenient to write it in a different way, using the following theorem. 

The absolute value of the flux of a uniform vector field over a surface, S, such that 
each field line intersects the surface at most once, is given by the product of the field 
intensity by the area of the projection, S pr , of the surface S on a plane orthogonal to 
the field. 

The proof follows from the definition of the flux, or Gauss's theorem. Using this result, 
the water mass can be written as 
L 

m{u) = -\}(u)\S pr (u). (7) 
u 

The time L/u clearly decreases monotonically with u. The projected wet surface 
S pr depends on the shape of the body and on its orientation with respect to the rain and 
the direction of motion — we cannot say anything in general about its dependence on u. We 
note also that the projection of the wet surface exactly overlaps with the projection of the 
whole body surface. 

The absolute value of j (see equation (2)) increases monotonically with u if v x < 0. 
We know, when travelling in the rain, that if we speed up, the intensity of the rain seems to 
increase. On the contrary, if v x > 0, the apparent rain density has a minimum for u = v x . 

The ratio |j \/(u) does not depend on the body: it is therefore connected with some general 
features of the problem and it is worthwhile to study its dependence on u. Let us write it in a 
more explicit form: 



|j| pJ(Vx ~ U) 2 + V 2 + V\ 



- = N\-) ~ 2 +1 =n(-) -2-COS0O + 1- (8) 
u u V \u/ u V \u/ u 

It is remarkable that this function, while monotonically decreasing for v x < 0, always has 

a minimum when v x > 0, that is for 

2 

V V 

« = - = —-. (9) 

V x COS #0 

A plot of |j \/(u) versus u/v for some values of #o is shown in figure 2. We can see that for 
= 30° the function has a well pronounced minimum, while for 60° the minimum is already 
hard to distinguish. 

This approach automatically solves the problem for all situations where the projected wet 
surface does not depend on u. Of course, the simplest example one can consider is a sphere, 
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1 2 3 4 5 6 

u/v 

Figure 2. Plot of the ratio |j|/(pw) (arbitrary units) versus u/v for the following values of the angle 
<9 between u and j : 30° and 60° (tailwind), 110° and 140° (headwind). 

for which the projected wet surface is always tvR 2 , R being the sphere's radius. But this is not 
the only case: as we have seen, when u changes, j rotates in a plane perpendicular to the £ 
axis, so every body with a symmetry axis along this direction is in the same situation. 

In conclusion, we can say that in all the cases where the projected wet surface does not 
depend on u. 

• An optimal speed exists, and its value is always given by equation (9), irrespective of the 
body shape (cylinder, cone, ellipsoid ...), under the only condition that the wind comes 
from behind. We note that the u value (9) is always greater than v x . Thus the optimal speed, 
when it exists, is not always equal to v x , as found by previous authors. 

• The optimal speed depends on v, and then on the drop size (since v obviously depends on 
v z ). 

• An optimal speed exists irrespective of the intensity of the wind, provided that it has a tail 
component. 

Of course, in general S pr does depend on u, and its effect on m(u) may be very different 
from one case to another. We can say that it is more relevant for bodies with a high degree of 
asymmetry with respect to rotations around the f axis. We can be more precise only for some 
particular, simple situations. 

5. A parallelepiped and a plane surface 

As said in the introduction, previous papers (with the exception of Ehrmann and Blachowicz 
[2], who considered a vertical cylinder) have so far considered a parallelepiped with edges 
parallel to Cartesian axes (we shall call it a 'vertical' parallelepiped) and a plane surface. Let 
us now direct our attention to these two cases. 

5.1. A vertical parallelepiped 

The case of a vertical parallelepiped has been discussed by several authors, but we shall restate 
the results briefly here with our notation for the sake of convenience. 
The rain flux results: 



<D(k) = p(S x \v x - u\ + Sy\v y \ + S z \v z \). 



(10) 
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Here S x denotes the area of each face perpendicular to x axis, and so on. As a consequence, 
the mass of the intercepted water is: 

m(u) = — (S x \v x - u\ + Sy\v y \ + S z |v z |). (11) 
u 

Depending on the value of v x , we must consider the following cases. 

• v x ^ 0: headwind or no wind. 

In this case \ v x — u\ = \v x \ + u. Then: 

m(u) = pL\S x + I. (12) 

This function decreases monotonically, so the best strategy is to move as quickly as 
possible. 

• v x > 0: tailwind. 

At 'low speed' (u < v x ), \v x — u\ = v x — u, so that 

m(w) = pLl -S x H ^ I . (13) 

This function, too, decreases monotonically, so it is never convenient to move at a speed 
lower than v x . But what about a speed higher than v x l 
At 'high speed' (u ^ v x ), \ v x — u\ = u — v x , then: 

m(u) = pLlS x + y —± I (14) 

and we must further differentiate between two subcases, depending on the sign of the 
expression: 

-S x \v x \+S y \v y \+S z \v z \. (15) 

When this expression is strictly positive, m(u) also decreases for u ^ v x , so it is convenient, 
once more, to move at a maximum speed. But when, on the contrary, it is negative, m(u) 
increases for u ^ v X9 an optimal speed does exist and is always: 

u opt = v x . (16) 
We can observe that in this case the optimal speed depends only on the wind speed and 
not on drop size. We note also that, when expression (15) vanishes, the mass of water 
intercepted by the body is equal to the water found in the volume LS X and is independent 
of u, provided only that u > v x . 
It is worthwhile to write down the condition that (15) be negative in two different ways: 
Sx I v x I + S y I Vy 

S x > (17a) 

Vx 

S X \V x \+Sy\Vy 

v x > (IV/?) 

^x 

Equation (17a) tells us that, for a given v, an optimal speed exists only for parallelepipeds 
with a large S x ; (17b) tells that, for a given parallelepiped, an optimal speed exists only if v x 
is large enough. 

5.2. A generalisation? 

For a parallelepiped, the rain flux is the sum of three terms, each relative to a component of j: 

*=\jx\S X +\jy\Sy+\j z \S z . (18) 

It is very tempting, at this point, to generalize this result to a body of any shape, by simply 
defining S x to be the projected area of the body on a plane normal to x axis, and so on, and 
this seems to be the line followed by previous papers. However, even though very appealing, 
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this idea is not correct, because the presence of the modulus in definition (5) does not permit 
splitting of the sum of the component products. 

We have already seen a simple example where equation (18) is not valid. In the case of a 
sphere, according to this formula the flux should be (\j x \ + \j y \ + \j z \)rcR 2 , while it turns out 
to be \i\7tR 2 . Below we shall see some more examples. 



5.3. A plane surface 

In order to develop a deeper understanding of the physics of the problem, let us consider a 
plane-oriented surface and A to be the vector area associated with it. We need to consider two 
more quantities: the angle, (p, between A and j, whose initial value (the angle between A and 
jo) will be denoted by cpo, and the angle, iff , between A and u. Without any loss of generality, 
we can always choose the orientation of A in such a way that A x is positive. Equation (6) then 
gives: 

v A 



L pL 
m(u) = -|j • A| = — |(v-u) • A| = pL 

u u 



-A, 



(19) 



From this relation we see the following. 

• When v • A = 0, that is when the surface is parallel to the rain, m does not depend on u, 
so the speed does not matter at all. The surface simply collects all the water contained in 
the volume LA X . In particular, when A x = 0, that is when the surface is in the 7t plane, it 
remains dry. 

• When and A x = 0, the way to minimize m is to let u be as high as possible. 

• Lastly, when v • A and A x are both non-zero, the only relevant parameter is the angle 
(po, which determines the sign of the dot product v • A. Then we have that when 
cpo is acute, it is always possible to choose a speed such that the surface does not 
get wet at all. This optimal speed is clearly 

v • A cosffo 

u op t = — — = v — • (20) 

A x cos y/o 

When, on the contrary, cpo is obtuse, we see from equation (19) that m(u) decreases 
monotonically, so the best strategy is again to move as quickly as possible (the surface 
should move in the opposite direction in order to avoid getting wet). 

It can be useful to plot expression (19), which is essentially the modulus of a vertically 
shifted hyperbola branch, for the case of an acute and an obtuse value of cpo (see figure 3). 

In conclusion, we can say that for a plane surface the choice of the best strategy depends 
uniquely, for a given rain velocity, on the orientation of the surface with respect to the rain 
(which determines the sign of the dot product v • A) and on the direction of motion (which 
determines if A x vanishes or not). 

Compared to the case of a vertical parallelepiped, for a plane surface the following 
differences can be highlighted. 

• When m(u) has a minimum, it is always zero, that is, the surface remains dry. 

• An optimal speed can also exist with a headwind. 

• There is no condition on v x \ an optimal speed can also exist with a light wind. 

• The optimal speed can be either lower or higher than v x . 

• The value of u opt also depends on the modulus of the rain velocity, and then on the drop 
size. 

The condition that the angle cpo is acute has a simple physical meaning: given our 
convention to orient A so that A x > (A 'points ahead'), it simply means that the rain 
wets the rear face of the surface. In this case, increasing u, the rotation of j will ultimately 
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\ 120° 




\ 60° 



0,5 1 1,5 2 2,5 3 

u/v 



Figure 3. Plot of m (arbitrary units) versus u/v for a plane surface in case of an acute (60°) and 
obtuse (120°) angle <po between v and A. 




Figure 4. A plane surface parallel to the y axis moves in the rain, in the absence of a lateral wind. 
(Left) We have a tailwind and the angle (po between the rain and A is obtuse (the rain wets the front 
face of the surface). Increasing u, j rotates clockwise and never becomes parallel to the surface, 
so m(u) does not have a minimum. (Right) We have a headwind and the angle (po is acute (the 
rain wets the rear face). When j rotates, there is a particular (optimal) speed for which it becomes 
parallel to the surface. 

lead to wetting of the front face, and so there will be a particular value of u for which the rain 
will be parallel to the surface. 

Now we can understand that the condition of having a tailwind is neither necessary nor 
sufficient to have a minimum in the function m(u). Figure 4 explains this point for the cases 
of v y = (no cross component of the wind) and A y = (surface parallel to the y axis). 

From the discussion of this point — which concerns a very simple geometrical situation — 
we can learn how complex the problem is and how we must pay attention to such aspects as 
the orientation with respect to the direction of motion and of the rain. 

6. From a plane surface to a solid body 

Of course, a solid body is not a plane surface, it has several surfaces with different orientations, 
so it is not surprising that results found in the two cases are different. For example, there is no 
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Figure 5. A parallelepiped with two opposite faces parallel to the xz plane. By analogy with the 
previously considered situation, we call ifr the angle between Ai and the x axis. We choose Aj and 
A3 so that their x component is positive. 



way for a solid body to remain dry in the rain. However, the 'transition' from a plane surface 
to a solid body, for physical reasons, must be smooth. 

Let us consider a parallelepiped with a generic orientation. First of all we observe that 
in a parallelepiped each pair of opposite faces behaves in the rain like one plane surface, so 
a parallelepiped, for our discussion, is equivalent to three mutually perpendicular rectangles. 
If we associate a vector area with each pair we get three mutually orthogonal vectors, that 
we shall call Ai, A2, and A3, and that we shall orient, as before, in such a way that their x 
component is positive. 

The equation that gives the water mass impinging on the parallelepiped is a simple 
generalisation of equation (19), which reduces to (11) for a vertical parallelepiped: 

/ v • Ai v • A2 

m(u) = pLi A ix + A 2x 

\ u u 



+ 



(21) 



This is the sum of three functions, each of which can have one of the two trends shown in 
figure 2. The study of this function is a very complex problem, so we shall focus our attention 
on a slightly simpler situation, that is an analogue of the 'bi-dimensional' case considered 
in figure 3: a parallelepiped with two faces perpendicular to the y axis (let us call A2 the vector 
area associated with this pair) and no cross component of the wind (v y = 0) (see figure 5). 

In this situation there is no flux on the faces perpendicular to the y axis. In order to shed 
some light on the transition from a plane surface to a solid body, we may plot function (21) 
for increasing values of the ratio A3/A1 (see figure 6). 

The plot has been obtained with the following parameters: #0 = 125° (headwind) and 
f = 50°. 

We can see how the thickness affects m(u) : it gradually makes the minimum shallower, and 
ultimately destroys it. The reason is clear; increasing the thickness has the effect of reducing 
the asymmetry of the projected surface with respect to rotations in the 71 plane (in this case, 
the xz plane). At the same time, it can be seen from this plot that there are cases where, even 
for a solid body, a minimum exists with a headwind, and where the optimal speed is lower 
than v x \ some of the differences between a plane surface and a parallelepiped, described in 
section 5.3, are just due to the orientation of the solid body. 
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Figure 6. Plot of m (au) versus u/\v x \ for a parallelepiped with ^ = 50° and rain falling at an 
angle 6q = 125° (headwind), for various values of the ratio A3 /A\ \ 0, 0.15, 0.30, 0.45. 




Figure 7. A vertical cylinder. The wet surface is highlighted. 



7. A cylinder 

7.1. Vertical cylinder 

Let us now consider a cylinder oriented vertically. In this case the wet surface S w is constituted 
by the upper base and by one half of the lateral surface. To be more precise, it is one of the two 
halves defined by a plane passing through the cylinder axis and perpendicular to the horizontal 
component of j, which we will denote j xy . The situation is shown in figure 7. 

On the upper base, the only contribution to the flux comes from j z . On the lateral surface, 
on the other hand, only j xy contributes. With regard to this latter term, it is evident that the 
associated projected surface S pr is a rectangle with base 2R and height H, R and H being, 
respectively, the radius and height of the cylinder. Then the rain flux is: 



d>0) = pi 7tR 2 \v z \ + 2RHJ(v x - u) 2 + v 2 



Consequently, the water mass is: 



m(u) = pLR 



(22) 



(23) 
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|v y |/v x = |v z |/v x = 4 H/R = 12 |v y |/v x = 0.5 |v z |/v x = 4 H/R = 12 
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u/v x u/v x 



Figure 8. Qualitative plots of m (arbitrary units) versus u/v x for a vertical cylinder for various sets 
of parameters. The upper left plot is identical to what can be obtained for a parallelepiped. 

This expression is functionally different from that obtained for a parallelepiped, showing 
one more case where equation (18) is not correct. Without performing a complete study of this 
function, we can obtain some information from it. 

First of all, we may note that for v x < this function is monotonically decreasing; a 
minimum may exist only with a tailwind. 

In case that v y = 0, the cylinder behaves just like a parallelepiped with A x = 2RH and 
A z = jtR 2 . In this particular case, then, equation (18) may be applied. A vertical cylinder was 
considered by Ehrmann and Blachowicz [2] who, however, did not take into account a cross 
component of the wind and so missed the difference with a parallelepiped. 

On the other hand, if v y ^ 0, we can plot the function (23) for various sets of parameters 
in order to get some idea of its general behaviour (see figure 8). 

Looking at these plots, we could conjecture the following. 

• Increasing the ratios \ v y \/v x or \ v z \/v x has the effect of making the minimum smoother and 
of shifting the optimal speed towards higher values. If these ratios are too high, a minimum 
may not exist at all. 

• The same effects can be obtained by decreasing the ratio H/R. If this ratio is too low, a 
minimum may not exist at all. 

8. Conclusion 

8.1. Results 

Starting from very general assumptions, an equation that describes how to calculate the mass 
of the water received by the body has been derived (we wrote it in two different ways: 
equations (6) and (7)). We have solved or studied this equation for a plane surface and for 
bodies with a simple shape. For a plane surface, we have found the following. 
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• An optimal speed exists subject to the sole condition that the rain wets the rear face of the 
surface, irrespective of the intensity of the wind and of the sign of v x . 

• The optimal speed is u opt = v • A/A*, so it depends on |v|, and then on the drop size. 

• When moving at this optimal speed, the surface does not get wet. 

• The optimal speed can have any value. 

For solid bodies, we have solved the equation in all cases where the area of the wet surface 
does not depend on u, that is when the body has a symmetry axis and this is perpendicular to 
the plane identified by v and u. In these cases, an optimal speed exists, provided that the wind 
has a component from behind; its value is u opt = v/ cosOq, where Oq is the angle between v 
and u, so in this case it also depends on the drop size. 

The case of a vertical parallelepiped had been already studied and solved by previous 
authors. We have extended this study by considering a slanted parallelepiped, and have found 
that in some cases it is possible that an optimal speed exists even with a headwind, and that its 
value is not always v x . 

Lastly, we have considered a vertical cylinder, showing that its behaviour is different from 
that of a vertical parallelepiped if there is a cross component of the wind velocity. 

Our study of the behaviour of solid bodies in the rain has revealed a wide range of 
situations, and general rules cannot be found. We can say that the presence of a tailwind seems 
to be a favourable, but not always a necessary condition for the existence of an optimal speed. 
In some cases, we have found that the value of u opt depends on the drop size, and in other 
cases it does not. 

8.2. Didactic considerations 

While the problem seems to be too complex to be proposed as a didactic exercise, some 
parts of it (e.g. the vertical parallelepiped, or the plane surface) could be suitable for didactic 
purposes at the undergraduate level. First of all, the problem is familiar to all, and as such 
stimulates students to apply themselves to find the solution, and also to gaze around with a 
'physical mind', looking for all the physical phenomena present in everyday life. Moreover, in 
my opinion it is instructive to see how concepts (like flux, current density, Gauss's theorem) 
turn out to be useful in other contexts. At the same time, it is beneficial to think about the 
fact that concepts are tools, and that they can and should be modified, if necessary, in order to 
adapt them to our purposes, as we have with the flux concept. 
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